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Field-effect  transistors  (FETs)  in  conventional  electronic  circuits  are  in  contact  with 
the  high-thermal-conductivity  substrate.  In  contrast.  FETs  in  novel  silicon-on-insulator 
( SOI)  circuits  are  separated  from  the  substrate  by  a  thermally  resistive  silicon-dioxide 
layer.  The  layer  improves  the  electrical  performance  of  SOI  circuits.  But  it  impedes 
conduction  cooling  of  transistors  and  interconnects,  degrading  circuit  reliability.  This 
work  develops  a  technique  for  measuring  the  channel  temperature  of  SOI  FETs.  Data 
agree  well  with  the  predictions  of  an  analytical  thermal  model.  The  channel  and  inter¬ 
connect  temperatures  depend  strongly  on  the  device  and  silicon-dioxide  layer  thick¬ 
nesses  and  the  channel-interconnect  separation.  This  research  facilitates  the  thermal 
design  of  SOI  FETs  to  improve  circuit  figures  of  merit,  e.g.,  the  median  lime  to  failure 
( MTF)  of  FET- interconnect  contacts. 


1  Introduction 

The  performance  and  reliability  of  electronic  circuits  arc  af¬ 
fected  by  temperature  fields  in  transistors  and  interconnects.  As 
transistor  dimensions  decrease,  thermal  conduction  within  a  few 
micrometers  of  these  heat  sources  governs  an  increasing  fraction 
of  the  transistor-to-coolant  temperature  difference  in  an  elec¬ 
tronic  system.  This  is  very  important  in  novel  silicon-on-insulator 
(SOI)  electronic  circuits,  where  transistors  are  separated  from 
the  substrate  by  a  thermally  resistive  silicon-dioxide  layer,  often 
fabricated  by  implanting  oxygen  ions  into  a  single-crystal  silicon 
wafer.  Figure  I  is  a  cross  section  of  a  SOI  field-effect  transistor 
( FET),  whose  common  dimensions  are  given  in  Table  1 .  Almost 
all  of  the  device  power  is  dissipated  in  the  channel.  The  electri¬ 
cally  insulating  implanted  layer  prevents  latchup  between  de¬ 
vices  and  reduces  the  parasitic  capacitance  of  the  transistor  due 
to  the  substrate,  facilitating  faster  circuits  (Colinge,  1991 ).  But 
the  implanted  silicon-dioxide  layer  has  a  low  thermal  conductiv¬ 
ity,  and  impedes  conduction  cooling  of  the  channel  through  the 
substrate.  The  SOI  circuit  designer  must  know  the  resulting  tem¬ 
perature  rise  in  interconnects  and  devices. 

McDaid  et  al.  (1989)  showed  that  the  temperature  rise  de¬ 
creased  the  drain  current  of  a  SOI  FET  for  given  gate  and  drain 
voltages  by  assuming  one-dimensional  heat  conduction  through 
the  implanted  layer  from  an  isothermal  FET.  Goodson  and  Flik 
( 1992)  predicted  the  temperature  field  in  a  SOI  FET  by  treating 
the  source,  drain,  gate,  and  interconnects  as  cooling  fins  for  the 
Joule-heated  channel.  They  showed  that  the  temperature  varies 
significantly  within  the  FET,  and  indicated  that  the  temperature 
rise  could  reduce  the  electromigration-limited  reliability  of  inter¬ 
connects.  For  application  of  SOI  circuits  below  77  K,  e.g.,  in 
hybrid  superconductor-semiconductor  circuits,  phonon-boundary 
scattering  was  shown  to  influence  the  temperature  field  strongly 
in  FETs. 

Experimental  confirmation  of  the  analysis  of  Goodson  and  Flik 
( 1992 )  requires  a  technique  for  measuring  temperature  locally  in 
transistors  with  spatial  resolution  comparable  to  the  channel 
length,  2Lk,  which  can  be  smaller  than  0.5  pm.  Lifka  and  Woer- 
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lee  (1990)  estimated  the  transistor-to-substrate  thermal  resis¬ 
tance  for  SOI  circuits  through  the  local  melting  of  a  coaling,  an 
approach  that  lacks  the  needed  spatial  resolution.  Bunyan  et  al. 
( 1992)  used  noise  thermometry  to  measure  temperatures  in  SOI 
transistors.  This  approach  requires  experimental  structures  very 
different  from  those  in  a  circuit,  yielding  an  impact  on  the  tem¬ 
perature  field  that  needs  to  be  assessed.  For  non-SOI  semicon¬ 
ductor  devices,  Brugger  (1991)  and  Ostermeier  et  al.  (1992) 
performed  temperature  measurements  with  submicrometer  spa¬ 
tial  resolution  using  micro-Raman  spectroscopy.  But  the  impact 
of  the  incident  radiation  on  the  performance  of  transistors  is  not 
known.  The  promising  method  of  Majumdar  et  al.  (1993)  and 
Lai  et  al.  ( 1993)  uses  an  atomic  force  microscope  to  determine 
local  temperature  fields  in  circuits.  Future  research  will  almost 
certainly  make  this  method  effective  for  transistors  with  channel 
lengths  much  less  than  I  qm. 

This  work  develops  a  technique  to  measure  the  channel  tem¬ 
perature  of  SOI  FETs  with  a  resolution  in  the  direction  of  current 
flow  of  0.32  pm.  the  FET  channel  length.  The  gate  serves  as  an 
electrical-resistance  thermometer  for  the  channel  temperature. 
This  approach  was  used  by  Mautry  and  Trager  ( 1990)  for  bulk 
( non-SOI )  circuits,  in  which  the  characteristic  length  scale  of  the 
temperature  field  is  much  larger  than  FET  dimensions.  The  gate 
is  not  isothermal  in  SOI  circuits,  in  which  the  high  thermal  re¬ 
sistance  of  the  implanted  layer  causes  a  large  fraction  of  the  chan- 
nel-to-coolant  temperature  difference  to  occur  within  transistors. 
In  the  present  work,  thermal  analysis  yields  the  average  channel 
temperature  from  the  electrical  resistance  measured  along  the 
nonisothermal  gate.  Room-temperature  data  are  compared  with 
predictions  of  the  analysis  of  Goodson  and  Flik  ( 1992).  Analysis 
estimates  the  impact  of  the  channel -interconnect  separation  and 
the  implanted-silicon-dioxide  layer  thicknesses  on  the  reliability 
of  FET -interconnect  contacts. 

This  work  helps  to  determine  the  effect  of  the  implanled-layer 
thermal  resistance  on  the  practical  potential  of  SOI  technology. 
The  experimental  method  developed  here  investigates  thermal 
conduction  processes  within  and  very  near  transistors  and  inter¬ 
connects.  These  conduction  processes  must  be  understood  before 
transistor  design  based  on  thermal  analysis,  i.e.,  device-level 
thermal  design,  can  accompany  traditional  electrical  design  of 
devices  to  yield  circuits  of  optimal  performance  and  reliability. 

2  Thermal  Analysis 

SOI  circuits  are  made  of  materials  with  very  different  thermal 
conductivities,  shown  in  Table  2.  Goodson  and  Flik  ( 1992)  used 
the  difference  in  the  conductivities  to  develop  a  simple  thermal 
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model  lor  thermal  conduction  in  SOI  circuits.  Section  2.1  reviews 
the  analysis  and  Section  2.2  applies  it  to  the  experimental  struc¬ 
ture. 

2.1  Steady-State  SOI  FET  Thermal  Model.  The  model  of 
Goodson  and  Flik  (1992)  is  for  steady-state  FET  operation, 
which  is  the  case  in  most  measurements  of  FET  electrical  prop¬ 
erties.  The  time  required  for  steady  state  to  be  achieved  is  near 
(</„):/(k„/C.,),  where  k„/C„  is  the  thermal  diffusivity  of  silicon 
dioxide,  yielding  about  200  ns  at  300  K.  The  model  is  also  useful 
for  the  case  of  steady-periodic  power  dissipation  in  clock-driven 
circuits  at  points  separated  from  the  channel  heater  by  at  least 
one  thermal  penetration  depth.  The  thermal  penetration  depth  in 
the  silicon  source  and  drain  is  approximately  (TkJC*)'n  =  0.4 
(itn,  where  k,,IC,,  is  the  thermal  diffusivity  of  heavily  doped  sil¬ 
icon,  and  r  =  5  ns  is  a  typical  clock  period. 

Heat  (low  from  the  tops  of  devices  and  interconnects  is  shown 
in  Section  2.2  to  be  negligible.  The  channel  was  modeled  as  an 
isothermal  heating  source.  This  neglects  the  complex  distribution 
of  heating  intensity  in  the  channel,  yielding  a  small  error  in  the 
average  channel  temperature  that  is  estimated  at  the  end  of  Sec¬ 
tion  2.1,  Variations  in  the  temperature  of  the  substrate -silicon 
dioxide  interface  are  small  compared  to  the  channel-temperature 
rise.  This  interface  was  assumed  to  have  the  uniform  tempera¬ 
ture  /'„. 

The  source,  drain,  gate,  and  interconnects  were  modeled  by 
Goodson  and  Flik  ( 1992)  as  one-dimensional  cooling  fins  for  the 


channel,  arranged  as  shown  in  Fig.  2.  The  heat  loss  per  unit  fin 
length  was  assumed  to  be  h, w, ( T,  -  T„),  where  7)  is  the  local 
fin  temperature  and  w,  is  the  fin  width’in  the  direction  normal  to 
heat  flow.  The  heat  transfer  coefficient  is  h,  =  'i'k.Jd,,,,  where  d,„ 
is  the  thickness  and  kjd,„  is  the  inverse  volume  resistance  of  the 
silicon-dioxide  layer  beneath.  The  dimensionless  function  'l'  was 
derived  by  Goodson  and  Flik  (1992)  to  account  for  two-dimen¬ 
sional  conduction  in  the  silicon  dioxide.  The  function  depends 


Nomenclature 


A  =  area,  m! 

A,  =  fin  cross-sectional  area,  m2 
a  =  characteristic  length  of  test 
structure,  m 

C  =  specific  heat  at  constant  vol¬ 
ume  per  unit  volume, 

J  m  '  K-‘ 

c,,  =  specific  heat  at  constant  pres¬ 
sure  per  unit  mass,  J  kg'1  K  1 
d  =  layer  thickness,  m 
dr  =  total  thickness  of  thermally 
grown  and  CVD  silicon  diox¬ 
ide,  m 

dK„  =  thickness  of  silicon-dioxide 
layer  between  gate  and  chan¬ 
nel,  m 

d„  =  thickness  of  implanted-silicon- 
dioxide  layer,  m 

di„  =  thickness  of  silicon  dioxide  be¬ 
tween  fin  and  substrate,  m 
£,  =  clectromigration  activation  en¬ 
ergy,  J 

Fr,(Y)  =  shape  function  for  temperature 
distribution  in  gate 
G  =  channel-to-air  thermal  conduc¬ 
tance,  W  K  '■ 

g  =  acceleration  due  to  gravity, 
m  s~2 

h  =  heat  transfer  coefficient,  W 
m  3  K  1 

ln  =  drain  current,  A 
/,,  =  gate  current,  A 
J  =  current  density,  A  m'2 
K,  =  electromigration  constant,  s 
k  =  thermal  conductivity, 

Wm  1  K  1 

k/,  =  Boltzmann  constant  = 

1.38  x  10  J  K  1 


k,  =  channel  thermal  conductivity, 
Wm  '  K  ' 

k,.  =  thermal  conductivity  of  silicon 
dioxide,  W  m"1  K  1 
k,  =  thermal  conductivity  of  lightly 
doped  silicon,  W  m  1  K  1 
L,,  =  separation  between  gate  and 
metal  interconnect,  m 
L,  =  half-length  of  gate  in  X  direc¬ 
tion,  m 

L,„  =  half-length  of  interconnect  be¬ 
tween  devices  in  X  direction, 
m 

MTF  =  median  lime  to  failure,  s 
m  =  ( hlkd)'n  =  inverse  thermal 
healing  length  of  fin.  m  ' 

P  =  device  power,  W 
Pmf  =  time-averaged  device  power, 

W 

R.  =  channel-to-substrate  thermal 
resistance  =  l/C.  K  W  1 
R,:  =  electrical  resistance  of  gale  in 
Y  direction,  Q 

R<£T)  =  gatc-clectrical-resistance  cali¬ 
bration  function,  $1 
T  =  temperature,  K 
T,  =  average  channel  temperature, 

K 

T(AY)  =  gate  temperature  distribution, 

K 

AT,,  =  average  temperature  rise  in 
gate.  K 

Tlll  llu,,  =  maximum  interconnect  temper¬ 
ature,  K 

T0  =  substrate  temperature,  K 
V,,  =  voltage  drop  along  gate,  V 
Vas  =  gate-source  voltage  drop,  V 
=  drain-source  voltage  drop,  V 
v  =  air  velocity,  m  s  1 


w,  =  fin  width  in  X-Y  plane, 
m 

w„  =  channel  width  in  Y  di¬ 
rection,  m 

w,  =  separation  between 
channel  and  gate  con¬ 
tacts  in  K  direction,  m 
w,„  =  width  of  metal  intercon¬ 
nect  in  Y  direction,  m 
X  =  coordinate  in  plane  of 
substrate,  m 

Y  =  coordinate  in  plane  of 
substrate,  m 

Z  =  coordinate  normal  to 
plane  of  substrate,  m 
Z|,  Z2,  Z\,  Z.,  =  constants,  Eqs.  ( I )— (4) 
0  =  coefficient  of  thermal 
expansion,  K  1 
c  =  emissivity 
H  =  viscosity,  kg  m  1  s  1 
p  =  mass  density,  kg  m  1 
a  =  Stcfan-Boltzmann  con¬ 
stant  =  5.67  x  10  K 
W  nT2  K  4 
t  =  clock  period,  s 

=  two-dimensional  con¬ 
duction  function 

Subscripts 

d  =  property  or  dimension 
of  source  and  drain 
/  =  property  or  dimension 
of  fin 

g  =  property  or  dimension 
of  gate 

m  =  property  or  dimension 
of  metal  interconnect 
o  =  property  or  dimension 
of  silicon  dioxide  layer 
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Table  1  Common  dimensions  of  SOI  FET  devices  and  those  of  the  test 
structures  used  here 


Dimensions 

In  FET 
(pm) 

In  Test 

Structure  (pm) 

0.4 

0.293  -  0.503 

additional -SiQj  thickness,  d. 

0.60 

interconnect  thickness,  dm 

mmi 

i 

device  thickness,  d,i 

0.08 

0.041  -0.177 

Kate  thickness.  </„ 

0.30 

0.29 

0.0055 

0.0055 

channel-interconnect  separation,  L,i 

0.5 

0.8  -  3.8 

gate  half-width,  Lr 

0.25 

0.16 

device  width,  wj 

0.8 

10 

channel-Kale  contact  separation,  w. 

2 

2 

only  on  the  ratio  of  the  width  of  the  fin  and  the  thickness  of  the 
thermally  resistive  layer  below  the  fin.  The  values  of  'I'  used  here 
range  between  1.016  and  1.028  for  conduction  from  the  source 
and  drain,  between  1.051  and  1.064  for  conduction  from  the  in¬ 
terconnects,  and  between  2.349  and  2.925  for  conduction  from 
the  gate.  The  largest  Biot  numbers  for  the  test  structures  in  the 
present  manuscript,  dfhftk „  where  k,  and  d,  arc  the  fin  conductiv¬ 
ity  and  thickness,  are  0.028  for  the  source  and  drain  fins,  0.029 
for  the  gate  fin,  and  0.00782  for  the  interconnect  fins.  The  fin 
thermal  healing  length  is  1/m  =  ( kAtlhw ,)|,J,  where  k  and  A, 
arc  the  thermal  conductivity  and  cross-sectional  area  of  the  fin. 
The  distance  from  a  heating  source  over  which  the  fin  tempera¬ 
ture  recovers  to  the  substrate  temperature  is  of  the  order  of  the 
healing  length.  The  healing  lengths  are  I  tm„  =  {kmd„Jh,„)'n  ~ 
7  \x m  for  the  interconnects,  Mm,,  =  ( k,,d,,lh,, )"1  ~  1  /rm  for  the 
source  and  drain,  and  l/m„  =  (kndt/hx),a  ~  1.5  /on  for  the  gate. 

The  length  in  the  X  direction  of  the  contact  between  the  in¬ 
terconnect  and  the  drain  is  approximately  3  /2m.  Because  this 
length  is  substantially  smaller  than  the  healing  length  of  the 
interconnect  fins,  Goodson  and  Flik  ( 1992)  neglected  the  dif¬ 
ferent  heat  transfer  coefficient  between  the  interconnect  and  the 
substrate  over  this  length.  A  single  fin  is  used  for  the  intercon¬ 
nect  and  the  contact  region  with  the  heat  transfer  coefficient 
that  prevails  when  the  interconnect  is  above  the  thermal  and 
CVD  silicon-dioxide  layers.  This  slightly  overestimates  the 
thermal  resistance  between  the  channel  and  substrate,  but  sub¬ 
stantially  simplifies  the  analysis. 

The  width  of  the  interconnect,  source,  and  drain  fins  in  the  Y 
direction,  10  /im,  is  larger  than  the  healing  length  in  these  fins. 
This  makes  possible  significant  temperature  variation  in  this 
direction,  which  is  not  considered  by  the  present  fin  analysis. 
Temperature  variation  due  to  conduction  from  the  sides  of  these 
fins  is  small  because  the  thermal  resistance  for  conduction  by  a 
unit  area  through  one  healing  length  in  the  fin,  approximately 
(d,„d,)'nl(k„k,)'n,  is  much  smaller  for  each  fin  than  the  ther¬ 
mal  resistance  for  conduction  from  a  unit  area  of  the  fin  side  to 
the  substrate,  approximately  d,„lk„.  But  temperature  variation 
in  the  Y  direction  due  to  conduction  through  the  gate  fin  is 
significant.  The  error  is  estimated  using  an  analytical  solution 
to  the  three-dimensional  thermal-conduction  equation  in  the  im¬ 
planted  silicon-dioxide  layer  and  the  channel,  source,  and  drain. 
The  difference  between  the  conductivity  of  the  channel  and  that 
of  the  source  and  drain  is  neglected  and  the  source  and  drain 
fins  are  assumed  to  be  very  long  in  the  X  direction.  Heal  gen¬ 
eration  in  the  channel  is  modeled  using  a  uniform  heat-flux 
boundary  condition  at  the  top  surface  of  the  channel.  Heat  loss 
to  the  gate  is  modeled  using  a  heat  flux  condition  from  the  top 
of  the  channel  near  the  edge,  w,,  -  d,,  s  |  Y  |  s  w,,.  The  average 
channel  temperature  rise  is  slightly  larger  than  those  calculated 
by  neglecting  temperature  variation  in  the  Y  direction  in  the 


Table  2  Thermal  conductivities  of  SOI  circuit  materials.  These  values  are 
discussed  in  Section  2.2.  The  channel  thermal  conductivity  was  not  used 
by  Goodson  and  Flik  (1992).  but  is  needed  in  Section  3.2  for  the  analysis 
of  the  experimental  data. 


Region  or 
Component 

Material 

Thermal 
conductivity  at 

T0  =  303  K 
(Wm-i  K->) 

substrate 

single-crystal  silicon, 

3  x  10‘ 5  boron  atoms  cm3 

ks=  148’ 

channel 

single-crystal  silicon 

6  x  101 7  boron  atoms  cm-3 

kc=  148" 

source  and 
drain 

single-crystal  silicon, 

1  x  1020  arsenic  atoms  cm*3 

V) 

NO 

II 

gate 

polysilicon, 

I  x  I020  arsenic  atoms  cm  3 

**  =  30” 

interconnect 

aluminum, 

I  mass-percent  silicon 

km  -  239* 

SOI  implanted 
layer 

silicon  dioxide, 
implanted 

*„=  1.40*** 

other 

insulating  layers 

silicon  dioxide, 
thermally  grown  and  CVD 

ka=  1.40"* 

‘Touioukianetal.  (1970)  ’*Tai  et  al.  (1988)  "'Sugawara  (1969) 


channel,  which  is  the  approximation  of  the  simple  multifin  anal¬ 
ysis  used  in  this  manuscript.  The  relative  error  in  the  predictions 
given  here  is  less  than  3  percent  for  a  channel  width  w,,  =  0.3 
jim  and  less  than  7  percent  for  w,,  =  10  /ini. 

The  devices  were  assumed  to  be  in  an  infinite  linear  array, 
each  connected  by  an  interconnect  of  length  2 Lm ,  and  each  dis¬ 
sipating  the  same  power  P.  This  idealization  resulted  in  an  esti¬ 
mate  of  the  worst-case  temperature  distribution  in  a  real  circuit 
for  a  given  value  of  the  device  separation,  2 Lm.  It  yielded  the 
two  planes  of  symmetry  shown,  which  were  adiabatic  bounda¬ 
ries.  The  temperature  and  location  in  the  interconnect  are  given 
by  T„,  and  the  parameter  ,v,„,  in  the  drain  by  T,,  and  x,h  and  in  the 
gate  by  T,  and  xs.  The  parameters  x„,,  x,h  and  xK  arc  not  related 
to  the  coordinates  X,  Y.  and  Z.  The  channel  temperature  is  T, . 
The  gate  is  separated  from  the  channel  by  a  silicon-dioxide  layer 
of  thickness  dc„  =  5.5  nm,  whose  thermal  resistance  is  negligible, 
yielding  Tt(xx  =  0)  =  T...  The  fin  equations  were  solved  by 
requiring  temperature  continuity  and  energy  conservation  at  the 
fin  interfaces,  yielding 

Tm  -  Tu  =  Z,  cosh  |m,„(L«  -  x,„)|  ( I ) 

T,i  -  T„  =  Z,  exp  [m, +  Z,  exp|  (2) 


1992) 
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Tk  -  r„  =  Z,  exp[-mc.ctl 


(3) 

The  variables  Z, ,  Z2 ,  Z, ,  and  Zi  were  determined  by  solving  a 
set  of  four  simultaneous  algebraic  equations,  given  in  matrix 
form  by  Goodson  and  Flik  ( 1992).  The  channel  temperature  is 
Tr  =  TK(xll  =  0)  =  Zi  +  T„.  The  maximum  temperature  in  the 
interconnect,  at  the  interconnect-device  contact,  is  T = 
T„Xx„,  -  0)  =  Z;  +  Z(  +  r„.  The  channel-to-substrate  thermal 
resistance  is  R,  =  Zi/P. 

The  isothermal-channel  approximation  does  not  account  for 
the  strong  spatial  dependence  of  the  rate  of  heat  generation  in  the 
channel  (e.g.,  Ostermcieret  al.,  1992;  Fushinobu  and  Majumdar, 
1993).  The  greatest  error  occurs  when  the  transistor  is  saturated, 
which  results  in  heat  generation  predominantly  near  the  channel- 
drain  interface.  An  upper  bound  for  the  error  is  estimated  using 
the  fin  method  of  Goodson  and  Flik  ( 1992)  and  the  assumption 
that  all  of  the  heat  generation  occurs  within  a  region  at  the  chan¬ 
nel-drain  interface  of  length  in  the  X  direction  two  orders  of 
magnitude  less  than  the  channel  length,  LJ 50.  Conduction  to  the 
gale  and  the  difference  between  the  channel  thermal  conductivity 
and  that  of  the  source  and  drain  arc  neglected.  This  approximate 
analysis  of  the  error  yields  a  highly  nonuniform  channel  temper¬ 
ature  rise.  But  the  average  channel-temperature  rise  differs  by 
less  than  0.5  percent  from  that  predicted  by  assuming  an  iso¬ 
thermal  channel.  The  analysis  of  Goodson  and  Flik  ( 1992)  is 
therefore  appropriate  for  predicting  the  average  channel  temper¬ 
ature  rise,  which  is  measured  by  the  experimental  technique  de¬ 
veloped  here.  In  contrast,  this  analysis  of  the  error  shows  that  the 
isothermal-channel  approximation  results  in  a  significant  under¬ 
prediction  of  the  temperature  rise  of  the  drain-interconnect  con¬ 
tact.  For  L,i  =  0.5  and  1 .5  /r m,  the  relative  errors  arc  1 3  and  1 1 
percent,  respectively. 

The  analysis  neglects  thermal  boundary  resistances.  Goodson 
et  al.  ( 1994)  measured  the  effective  thermal  conductivities  for 
conduction  normal  to  the  implanted  silicon-dioxide  layers  in  the 
present  study  bounded  below  by  silicon  and  above  by  aluminum. 
The  data  agree  closely  with  the  conductivity  of  bulk  silicon  di¬ 
oxide,  making  a  significant  thermal  boundary  resistance  between 
the  implanted  layer  and  the  silicon  dioxide  unlikely.  Data  arc 
needed  for  the  boundary  resistance  between  the  buried  silicon 
dioxide  and  the  silicon  device,  as  well  as  the  resistances  at  the 
boundaries  of  the  thermally  grown  silicon  dioxide.  Goodson  et 
al.  ( 1993)  measured  the  effective  conductivities  for  conduction 
normal  to  chcmical-vapor-deposited  layers  fabricated  using  the 
same  method  as  those  in  the  present  study.  They  obtained  upper 
bounds  for  the  thermal  boundary  resistance  between  the  alumi¬ 
num  and  the  silicon  dioxide  that  would  not  significantly  change 
the  values  of  h ,  for  the  interconnects  in  the  present  work.  The 
thermal  resistances  that  impede  conduction  from  the  channel  to 
the  gate  and  from  the  source  and  drain  to  the  interconnects  are 
unknown.  If  significant,  these  resistances  will  cause  the  analysis 
here  to  underpredict  the  channel-to-substrate  thermal  resistance. 

2.2  Application  to  the  Experimental  Test  Structure 

Thermal  Conductivities.  Table  2  gives  the  thermal  conduc¬ 
tivities  used  in  the  model.  Goodson  et  al.  ( 1994)  measured  the 
thermal  conductivity  of  implanted-silicon-dioxidc  layers  in  SOI 
wafers  near  room  temperature.  The  data  agreed  within  the  ex¬ 
perimental  error  with  the  value  recommended  for  bulk  amor¬ 
phous  silicon  dioxide  (Sugawara,  1969),  so  the  bulk  value  is 
used  here.  The  bulk  value  is  also  used  for  thermally  grown  silicon 
dioxide,  supported  by  the  data  of  Goodson  et  al.  ( 1993).  The 
thermal  conductivities  reported  for  chemical-vapor-deposited 
(CVD)  silicon-dioxide  layers  do  not  agree,  but  are  in  general 
less  than  the  bulk  value  ( Schafft  et  al.,  1989;  Brotzenetal.,  1992; 
Goodson  et  al.,  1993).  Due  to  the  lack  of  a  consensus  among  the 
data  for  CVD  silicon-dioxide  layers,  the  bulk  value  is  used  in 
this  case  as  well.  The  resulting  error  is  very  small  because  the 
thermal  resistance  of  the  implanted  silicon-dioxide  layer  is  much 
more  important  than  that  of  the  CVD  layers. 
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The  source  and  drain  are  single-crystal  silicon  doped  with  ap¬ 
proximately  1  x  10’"  arsenic  atoms  cm  '.  The  most  appropriate 
existing  data  are  for  bulk  single-crystal  silicon  doped  with  1 .7  X 
10;"  phosphorus  atoms  cm  1  (Touloukian  et  al.,  1970).  There 
are  no  thermal  conductivity  data  available  for  the  gates  in  the 
present  research,  which  were  polysilicon  heavily  doped  with  ar¬ 
senic  atoms.  Tai  et  al.  ( 1988)  and  Volklein  and  Baltes  ( 1992) 
measured  thermal  conductivities  near  30  W  m  1  K  1  in  polysil¬ 
icon  layers  heavily  doped  with  phosphorus.  The  data  arc  consis¬ 
tent  with  the  thermal  tliffusivity  data  of  Mastrangelo  and  Muller 
( 1988),  also  for  phosphorus-doped  polysilicon  layers.  The  con¬ 
ductivity  =  30  W  m  1  K'1  is  used  here.  The  conductivity 
above  300  K  of  silicon  with  less  than  10,K  dopant-atoms  cm  ' 
differs  little  from  that  of  intrinsic  silicon  (Touloukian  et  al., 
1 970 ) ,  which  is  used  for  k,  and  k,  here.  The  thermal  conductivity 
of  aluminum  layers  containing  I  mass  percent  of  silicon  has  not 
been  measured  directly.  The  thermal  conductivity  calculated  us¬ 
ing  the  Wiedemann-Franz  law  (Kittel,  1986)  and  the  electrical 
resistivity  measured  here  of  these  layers  is  within  4  percent  of 
the  thermal  conductivity  recommended  for  bulk  aluminum  (Tou¬ 
loukian  et  al.,  1970),  so  the  bulk  value  is  used  here. 

The  substrate  temperature  during  the  measurements  was  T„  = 
303  K,  and  the  largest  channel  temperature  was  T,  =  403  K.  The 
thermal  model  in  Section  2.1  neglects  the  temperature  depen¬ 
dence  of  the  thermal  conductivities  of  the  SOI  FET  materials. 
This  is  a  good  approximation  between  303  and  433  K  for  silicon 
dioxide  and  aluminum,  whose  bulk  thermal  conductivities  vary 
by  less  than  13  and  2  percent  in  this  range,  respectively.  But  the 
thermal  conductivity  of  the  heavily  doped  silicon  source  and 
drain  varies  more  significantly  in  this  temperature  range.  To  help 
overcome  this  difficulty,  the  temperature  (T),  +  T,)l 2,  averaged 
for  all  of  the  data,  is  used  when  interpolating  k,,. 

Dimensions.  Figure  3  compares  the  experimental  test  struc¬ 
ture  with  a  FET  device.  The  measured  dimensions  of  the  test 
structures  arc  given  in  Table  1 .  The  interconnect  lengths  are  very 
long  compared  to  1  lm,„ ,  so  that  L„,  =  is  used.  The  gate  of  the 
FET  in  the  test  structure  extends  out  from  both  sides  of  the  chan¬ 
nel  to  interconnect  contacts  in  the  Y  direction.  This  results  in  a 
plane  of  symmetry  normal  to  the  Y  axis,  i.e.,  an  X-Z  plane,  at  Y 
=  0.  The  temperature  field  in  each  half  of  the  test  structure 
is  predicted  by  the  thermal  model.  The  thermal  analysis  is  ap¬ 
plied  using  w,  =  (wJ)la-UBJClure/2,  *v,„  =  (»v,„)*,,M,„,u,,./2,  and  P  = 
(  P  )tCM  MniKlUU'f  2. 


a)  FET  DEVICE.  b)  TEST  STRUCTURE. 

Fig.  3  Top  views  of:  (a)  FET  device;  (b)  test  structure 
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Parameter  Uncertainties.  An  uncertainty  in  the  predictions 
of  the  thermal  analysis  results  from  the  use  of  parameters,  e.g., 
thermal  conductivities  and  dimensions,  which  may  differ  from 
those  in  the  test  structure.  This  uncertainty  is  estimated  using  the 
sum-of-squares  technique  (e.g.,  Holman,  1984)  and  the  analyti¬ 
cal  model.  The  largest  expected  error  is  ±10  percent.  The  un¬ 
certainties  in  k,i  and  d,,  are  the  largest  contributors. 

Heat  Transfer  to  Ambient  Air.  The  test  structure  is  exposed 
lo  ambient  air,  but  heat  transfer  to  the  air  is  neglected.  This  is 
justified  by  the  small  value  of  the  channel-to-air  thermal  conduc¬ 
tance  compared  to  the  channcl-to-subsirate  thermal  conductance, 
which  is  predicted  in  this  work  to  be  \/R,  ~  0.5  -  2  X  10  4 
W  K  1 .  The  channel-to-air  thermal  conductance  is  of  the  order 
of  that  from  an  isothermal  disk  of  radius  a  on  the  boundary  of  a 
semi-infinite  medium  of  conductivity  k.  G  =  4 ak  (Carslaw  and 
Jaeger,  1959).  Using  a  =  w,,/2  and  the  room-temperature  con¬ 
ductivity  of  air  yields  C  =  5.2  X  10  7  W  K"1.  An  order-of- 
magnitude  analysis  of  the  momentum  equation  estimates  the  air 
velocity  near  the  device  due  to  buoyancy  forces  ( Rohsenow  and 
Choi,  1961 ),  u  =  pga2P(T,  —  T„)/p,  where  g  is  the  acceleration 
due  to  gravity,  and  p  is  the  density,  p  is  the  viscosity,  and  0  ~ 
I /Tu  is  the  approximate  coefficient  of  thermal  expansion  of  the 
air.  Using  room-temperature  properties  and  T,  -  T»  =  130  K 
yields  u  =  6.6  pm  s'1.  The  thermal  conductance  contributed  by 
the  air  motion  is  of  the  order  of  G  =  vira2pc,,  =  6.2  X  10"" 
W  K  1 ,  where  c(,  is  the  specific  heat  per  unit  mass  at  constant 
pressure  of  air.  The  thermal  conductance  due  to  radiation  is  of 
the  order  of  G  =  4(oT»A ,  where  t  is  the  cmissivity  of  the  surface, 
cr  is  the  Stefan- Boltzmann  constant,  and  A  is  the  area  of  the 
emitting  surface.  Using  e  =  1  and  A  =  7r  a2  yields  G  = 
5.0  X  10  10  W  K  1  at  room  temperature. 

3  Channel-Temperature  Measurement  Technique 

This  section  describes  the  technique  for  measuring  the  channel 
temperature  of  SOI  FETs.  Section  3.1  describes  the  apparatus 
and  the  general  procedure,  and  Section  3.2  calculates  the  channel 
temperature  from  the  measured  gate  resistance.  Section  3.3  de¬ 
termines  the  experimental  uncertainty. 


TC(Y)  =  T„  +  AT,;F,.(Y) 

where  A T„  is  the  average  gate-temperature  rise  from  and 
F,,{Y)  is  a  shape  function  of  average  value  unity  that  is  defined 
for  |  f  |  <  wr  +  w,,l 2.  For  each  measured  R„,  A T„  is  determined 
iteratively  using 

Fe.  =  f  lRc(T0  +  A TcFt!(Y))LM„  .  "  (5) 

2  2  wr  +  w,, 

The  thermal  resistance  of  the  silicon-dioxide  layer  between  the 
channel  and  gate,  =  1.2  x  lO'KW’1,  is  small 

compared  to  the  thermal  resistance  for  conduction  along  the  gate 
to  the  contact  in  the  Y  direction,  wr/(2LxdJtK)  =  7.31  X  10s  K 
W  This  means  that  the  channel-  and  gate-temperature  distri¬ 
butions  are  almost  identical  for  |  Y  |  <  w,,/2.  The  average  channel 
temperature  is 

dY 

7;.=  r„+  ATaFa(Y)  —  (6) 

J-«w«  wi 

The  channel-to-substrate  thermal  resistance  is  R,  =  (T,  -  T»)/P 
=  (T.~  To)/a,>Vlu). 

Two  shape  functions  are  now  developed,  from  which  Eqs.  (5) 
and  ( 6 )  yield  upper  and  lower  bounds  for  T,  for  a  given  R„ .  Each 
shape  function  must  be  even,  due  to  the  test-structure  symmetry 
about  Y  =  0,  and  continuous.  Because  of  the  large  width  in  the 
X  direction  of  the  gate-interconnect  contacts,  4  //m,  compared  to 
the  gate  width,  2LS  =  0.32  //m,  the  contacts  are  very  nearly  iso¬ 
thermal  at  the  substrate  temperature,  7j,.  This  yields  the  boundary 
conditions  Fa(Y)  =  0  at  Y  =  ±(wr  +  w,,/2).  For  a  given  R„,  T, 
calculated  using  Eqs.  (5)  and  (6)  increases  with  the  difference 
between  unity  and  the  average  of  Ft,(Y)  over  -tv,,/2  <Y<  w,,l 
2,  i.e.,  with  the  assumed  difference  between  the  average  gate  and 
channel  temperatures.  Shape  I,B  assumes  a  linear  temperature 
profile  in  the  gate  segments  not  over  the  channel,  which  neglects 
conduction  down  through  the  buried  silicon-dioxide  layer,  and 
an  isothermal  channel.  Both  assumptions  underestimate  the  dif¬ 
ference  between  the  average  channel  and  average  gate  tempera¬ 
tures,  yielding  a  lower  bound  for  T, .  Shape  LB  is 


3.1  Apparatus  and  Procedure.  Figure  3  shows  the  exper¬ 
imental  structure.  The  electrical  resistance  of  the  gate  depends 
strongly  on  temperature.  It  serves  as  an  electrical-resistance 
thermometer.  The  calibration  measures  the  gate  electrical  resis¬ 
tance,  Rft ,  as  a  function  of  temperature  when  there  is  no  drain 
current,  i.e.,  when  the  gate  is  isothermal,  yielding 
[/ffiCDLntanu™.-  The  substrate  temperature  is  controlled  using 
a  Temptronic  Model  TP38B  chuck,  a  copper  disk  with  a  di¬ 
ameter  of  88.9  mm  and  a  thickness  of  19.1  mm,  to  which  the 
wafer  is  secured  by  suction.  A  thermocouple  with  one  junction 
soldered  to  the  chuck  surface  measures  the  chuck  temperature. 
The  chuck  is  maintained  at  the  temperature  Tn  and  the  gate 
resistance  is  measured  for  varying  values  of  the  drain-source 
voltage  drop,  V„s,  and  the  gate-source  voltage  drop,  Vcs.  i.e., 
for  several  different  device  powers,  P  = 

The  average  gate  temperature  is  defined  as  that  of  the  gate 
segment  whose  resistance  is  measured,  i.e.,  the  segment  between 
the  voltage  contacts.  The  average  channel  temperature  is  T,, 
which  is  shown  in  Section  3.2  to  be  very  well  approximated  by 
the  average  temperature  of  the  gate  segment  over  the  channel. 
The  average  gate  temperature  considers  the  gate  segments  not 
over  the  channel  heater,  and  is  less  than  T. .  The  FET  gate-tem¬ 
perature  variation  is  more  important  in  a  SOI  wafer,  where  most 
of  the  temperature  drop  occurs  within  a  few  micrometers  of  the 
channel  due  to  the  implanted  layer,  than  in  a  normal  substrate, 
where  the  temperature-drop  length  scale  is  the  thickness  of  the 
substrate,  i.e..  a  few  hundred  micrometers.  This  temperature  vari¬ 
ation  must  be  considered  when  calculating  T,  from  R„. 

3.2  Temperature  Distribution  in  the  Gate.  The  gate  tem¬ 
perature  variation  in  the  Y  direction  is 


w,,  +  wr 


F<;(Y) 


+  2 wr  [(w,  +  wj 2)  -  |F|| 
wri  +  w,  w. 


(7) 


y<m<»V  +  y  (8) 

An  overestimate  of  the  difference  between  the  average  gate 
and  channel  temperatures  requires  an  overestimate  of  the  tem¬ 
perature  drop  between  the  center  and  the  edge  of  the  channel  in 
the  Y  direction.  This  is  calculated  by  isolating  the  channel  and 
gate  from  the  source  and  drain,  in  which  thermal  conduction 
reduces  the  channel-temperature  variation  in  the  Y  direction.  The 
gate  and  channel  are  grouped  together  as  a  composite  fin,  which 
meets  a  fin  of  different  internal  properties  and  heat  transfer  co¬ 
efficient  at  T  =  w, ,12.  Solving  the  heat  equation  in  the  two  fins 
yields  shape  UB,  which  has  a  larger  average  in  the  channel  region 
than  shape  LB.  Values  of  T,  -  T„  calculated  using  shapes  UB 
and  LB  differ  by  less  than  8  percent,  and  the  simpler  shape  LB 
is  used  here,  Eqs.  (7)  and  (8).  This  function  does  not  describe 
the  temperature  distribution  in  the  channel.  Rather,  it  is  a  shape 
function  which,  when  used  in  Eqs.  (5)  and  (6),  yields  a  value 
for  T,  close  to  the  average  channel  temperature.  The  difference 
between  the  upper  and  lower  bounds  is  used  in  the  experimental- 
uncertainty  analysis. 

3.3  Experimental  Uncertainty.  The  uncertainty  in 
R.  =  (T,  -  T„)/P  has  three  significant,  independent  com¬ 
ponents:  (a)  There  is  a  relative  uncertainty  of  4  percent  in 
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Fig.  4  Channel-temperature  data  tor  SOI  FETs  with  varying  device  thick¬ 
nesses,  compared  with  data  for  a  bulk  (non-SOI)  FET 


T,  -  Tn  due  to  the  error  in  the  substrate-temperature  change 
measured  by  the  chuck  thermocouple  (Goodson  et  al., 
1994).  (b)  A  relative  uncertainty  of  6.6  percent  in  T.  -  T„ 
is  due  to  the  measurement  of  R(,.  (c)  A  relative  uncertainty 
of  8  percent  in  T,  -  Tu  is  due  to  the  approximate  shape 
function  for  the  temperature  profile  in  the  gate,  as  shown  in 
Section  3.2.  The  total  relative  uncertainty  in  /?.  is  ±  1 1  per¬ 
cent,  determined  using  the  sum-of-squares  technique  (Hol¬ 
man,  1984). 

4  Results  and  Discussion 

Section  4.1  compares  data  for  the  channel-to-substraic  thermal 
resistance  of  SOI  FETs  with  predictions  of  the  thermal  analysis 
of  Section  2.  The  thermal  analysis  is  used  in  Section  4.2  to  es¬ 
timate  the  influence  of  the  implanted  silicon-dioxide  layer  on  the 
reliability  of  highly  integrated  SOI  circuits. 

4.1  Channel-to-Substrate  Thermal  Resistance.  Channel- 
temperature  measurements  are  performed  on  SOI  test  structures 
with  varying  values  of  L,,,  d,,,  and  d„.  Test  structures  fabricated 
from  bulk  (non-SOI)  wafers  are  measured  for  comparison.  The 
device  voltages  satisfy  0  V  <  Vns  <  3  V  and  V6i  =  2  V  and  2.5 
V,  which  are  typical  operating  conditions.  The  device  powers 
vary  between  3  and  14  mW.  and  the  values  of  T,  -  T„  vary 
between  5  and  130  K.  The  power  dissipated  in  the  gate  electrical- 
resistance  thermometer  is  at  least  two  orders  of  magnitude 
smaller  than  the  device  power. 

The  channel-to-substrate  thermal  resistance,  R, ,  varies  by  less 
than  the  experimental  uncertainty  for  varying  powers  from  a  sin¬ 
gle  device,  as  shown  in  Fig.  4  for  three  SOI  devices  and  one  bulk 
device,  i.e.,  a  device  fabricated  in  a  conventional  wafer,  which 
lacks  the  implanted  silicon-dioxide  layer  in  SOI  wafers.  The  data 
for  each  device  for  varying  powers  fall  near  a  line  originating  at 
P  =  0  and  Tn  =  303  K,  whose  slope  is  R, .  Values  of  R,  for  the 
SOI  devices  are  as  much  as  10  times  larger  than  R,  for  the  bulk 
device,  due  to  the  thermal  resistance  of  the  implanted  silicon- 
dioxide  layer.  The  value  of  R,  decreases  with  increasing  SOI 
device  thickness,  dH.  In  Figs.  5  and  6,  each  data  point  is  the 
average  of  the  values  of  R,  measured  in  a  single  test  structure. 
Uncertainty  bars  are  only  given  for  selected  data  to  avoid  clut¬ 
tering  of  the  figures.  Each  data  point  is  the  average  of  the  values 
of  R.  measured  in  a  single  test  structure.  Uncertainty  bars  are 
only  given  for  selected  data  to  avoid  cluttering  of  the  figures. 

Figure  5  shows  that  the  sensitivity  of  R,  to  the  device  thickness 
is  predicted  by  the  thermal  model  of  Section  2.  Increasing  d,, 
reduces  the  channel  temperature  for  a  given  power.  The  param¬ 
eters  k,i  and  d,t  are  not  independent,  but  always  appear  as  a  prod¬ 
uct  in  the  solution  for  the  temperature  distribution  (Goodson  and 
Flik,  1992).  Thus,  the  channel  temperature  is  also  sensitive  to 
k,<,  which  depends  on  the  doping  level  in  the  source  and  drain. 
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Fig.  5  The  channel-to-substrate  thermal  resistance,  Re,  as  a  function  of 
the  device  thickness 


The  agreement  is  excellent  considering  the  uncertainties  of  the 
thermal  conductivities  and  dimensions  used  in  the  analysis.  Sec¬ 
tion  2.2  estimated  that  the  potential  relative  error  in  the  predic¬ 
tions  due  to  these  uncertainties  is  ±10  percent.  The  differences 
in  the  experimental  data  at  device  thicknesses  near  45  pm  in  Fig. 
5  do  not  indicate  a  poor  repeatability  of  the  experimental  method. 
Each  data  point  is  for  a  different  experimental  structure.  Because 
R,  depends  sensitively  on  the  structure  dimensions,  the  uncer¬ 
tainly  in  the  measurement  of  structure  dimensions  can  yield  data 
on  a  single  graph  that  are  apparently  inconsistent  when  only  the 
uncertainty  in  the  measurement  of  R,  is  considered. 

Figure  6  shows  the  dependence  of  R,  on  the  implantcd-silicon- 
dioxide  layer  thickness.  The  data  support  the  predictions  of  the 
model,  and  indicate  that  R,  is  as  sensitive  to  d,,  as  it  is  to  d„.  This 
is  in  contrast  to  the  predictions  of  McDaid  ct  al.  ( 1989),  whose 
model  assumed  that  R,  is  independent  of  d,,.  By  modeling  one¬ 
dimensional  conduction  in  the  implantcd-silicon-dioxidc  layer, 
these  researchers  predicted  that  R,  =  d„/(Ak„),  where  A  is  the 
device  area  in  the  X -Y  plane.  This  neglects  the  spreading  of  the 
temperature  profile  into  the  source  and  drain  fins  with  increasing 
d„,  and  is  not  consistent  with  the  data.  This  can  be  remedied  by 
a  simple  scaling  analysis.  The  area  in  the  source  and  drain  with 
significant  temperature  rise  is  of  the  order  of  A  =  2 wl(/m,,,  where 
the  thermal  healing  length  1/m,,  is  approximately  (k,id,id„IK.)'n ■ 
Using  R,  =  d„HAk„)  with  this  expression  for  A  yields 


which  is  in  qualitative  agreement  with  the  data  in  Fig.  6.  This 
shows  that  R,  is  roughly  proportional  to  d'„n ,  and  that  the  sen¬ 
sitivity  of  R,  to  d„  and  d,,  is  similar,  i.e.,  halving  d„  and  doubling 
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Fig.  6  The  channel-to-substrate  thermal  resistance  as  a  function  of  the 
implanted  silicon-dioxide  layer  thickness 
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cl, i  have  the  same  impact,  if  all  other  parameters  are  held  constant. 
Equation  (9)  is  valid  only  when  the  thermal  healing  length  in  the 
source  and  drain  is  smaller  than  the  channel -interconnect  sepa¬ 
ration,  i.c.,  I/m,,  <  L,,.  Otherwise,  there  is  significant  heat  con¬ 
duction  into  the  interconnects,  which  are  more  effective  fins  than 
the  source  and  drain  because  of  their  large  thickness  and  high 
thermal  conductivities. 

4.2  lilcctromigralion-Limited  Reliability  of  FET-Intercon- 
nect  Contacts.  The  FET-interconncct  contact  temperature  in¬ 
creases  as  L,i  is  reduced,  because  this  brings  the  contact  nearer 
to  the  channel  heat  source.  This  reduces  the  reliability  of  the 
contact,  whose  electromigration-limited  mean  time  to  failure 
(MTF)  decreases  with  increasing  temperature.  This  section  es¬ 
timates  the  reduction  of  the  MTF  of  FET-interconnect  contacts 
due  to  the  implanted  layer  in  SOI  circuits  and  demonstrates  the 
potential  for  overcoming  this  problem  through  transistor-level 
thermal  design.  Electromigration  is  the  motion  of  metal  atoms  in 
the  interconnect  in  the  direction  of  electron  flow  due  to  electron- 
lattice  momentum  transfer.  This  causes  interconnect  failure  due 
to  void  formation,  particularly  near  FET-interconnect  contacts, 
where  the  flux  of  metal  atoms  diverges.  The  temperature  depen¬ 
dence  of  the  MTF  limited  by  electromigration  is  (c.g..  Black, 
1967) 

MTF  =  /Ccxp(^)  (10) 

where  F„.  is  the  activation  energy  for  atomic  diffusion, 
k„  =  1.38  x  10  3'  J  K  1  is  the  Boltzmann  constant.  T  is  the 
temperature,  and  K,.  is  a  function  of  the  electrical  current  density, 
the  geometry,  and  the  microstructure  and  purity  of  the  metal. 
Equation  (10)  agrees  well  with  data  for  FET-interconncct  con¬ 
tacts  if  E,  =  0.5  eV  =  8  X  10  J  is  used  (Chem  cl  al.,  1986). 

The  FET  device  dimensions  in  Table  1  arc  used  with  an  in¬ 
terconnect  length  between  devices  of  2  L„,  =  4  pm.  The  FETs 
experience  steady-periodic  heating  with  pulses  of  1 .61  W  for  one 
tenth  of  each  clock  cycle,  where  1 .61  W  is  the  steady-state  power 
of  a  SOI  device  with  these  dimensions  (Woerlec  et  al.,  1989). 
The  analysis  in  Section  2  provides  a  good  estimate  of  the  nearly 
steady-state  FET-interconncct  contact  temperature  if  the  time- 
averaged  power  is  used,  /*«,„  =  0.161  W.  The  resulting  contact 
temperature  rise  is  less  than  8  K,  and  depends  strongly  on  d„  and 
L,,.  Figure  7  uses  Eq.  ( 10)  to  show  the  ratio  of  the  MTF  for  FET  - 
interconnect  contacts  in  a  SOI  circuit  to  that  for  contacts  in  a 
bulk  circuit.  These  predictions  assume  that  the  substrate  temper¬ 
ature  in  the  SOI  ease  is  equal  to  that  in  the  bulk  case.  The  dif¬ 
ference  between  the  two  contact  temperatures  is  due  to  the  ther¬ 
mal  resistance  of  the  implanted  silicon-dioxide  layer.  The  MTF 
increases  as  L,,  is  increased,  because  the  interconnect  moves  away 
from  the  channel  heater.  It  may  be  possible  to  improve  circuit 
reliability  by  increasing  L,,_  but  this  must  be  weighed  against  the 
need  for  compact  devices.  Reducing  d„  also  increases  the  MTF, 
because  this  reduces  the  contact  temperature. 

Figure  7  provides  the  type  of  information  needed  to  make  de¬ 
cisions  effectively  during  the  design  of  SOI  circuits,  but  the  pre¬ 
dictions  are  very  approximate.  Because  Eq.  ( 10)  has  only  been 
experimentally  verified  using  accelerated  testing,  i.e.,  the  use  of 
electrical  current  densities  and  temperatures  which  are  higher 
than  those  found  in  an  operating  circuit.  Fig.  7  can  at  best  show 
the  expected  trends.  In  a  real  device,  the  rate  of  heat  generation 
will  be  greatest  near  the  drain  contact,  causing  this  contact  to 
have  the  lowest  MTF.  But  the  isothermal-channel  approximation 
of  Goodson  and  Flik  ( 1992)  yields  identical  temperatures  in  the 
source  and  drain  contacts.  As  discussed  in  Section  2. 1 ,  this  un¬ 
derpredicts  the  temperature  rise  of  the  drain  contact  by  as  much 
as  13  percent.  The  resulting  error  in  the  normalized  MTF  in  Fig. 
7  is  less  than  10  percent. 

5  Conclusions 

The  analysis  of  Goodson  and  Flik  ( 1992)  lends  to  underpredict 
the  data  presented  here.  This  is  due  in  part  to  the  assumption  of 

580 /Vol.  117,  AUGUST  1995 


CHANNEI.-INTERCONNECT  SEPARATION, 
<nm) 


Fig.  7  Predicted  dependence  of  the  median  time  to  failure  (MTF)  of  FET- 
interconnect  contacts  on  the  channel-interconnect  separation 


an  isothermal  channel,  which  results  in  about  a  7  percent  under¬ 
prediction  of  the  channel-to-substrate  thermal  resistance,  as  dis¬ 
cussed  in  Section  2. 1 .  The  underprediction  may  be  due  in  part  to 
thermal  boundary  resistances  in  the  structure,  which,  as  discussed 
in  Section  2.1,  arc  not  considered  by  the  analysis.  Another  im¬ 
portant  possibility  is  that  the  thermal  conductivities  of  the  inter¬ 
connects  and  source  and  drain,  which  have  not  been  measured, 
may  be  less  than  the  bulk  values  used  here.  Better  agreement  will 
require  more  research  on  thermal  conductivities  and  boundary 
resistances  in  circuits. 

The  mean  time  to  failure  of  a  circuit  and  the  channel  mobility 
are  important  to  the  design  of  circuits.  They  are  affected  by  the 
channel  and  interconnect  temperatures,  which  are  shown  here  to 
depend  strongly  on  design  parameters,  e.g.,  *,,,  d,,,  d„,  and 
Some  of  these  parameters  also  affect  the  electrical  performance 
of  the  device  directly.  In  order  to  achieve  circuits  of  optimal 
performance  and  reliability,  design  for  electrical  performance 
should  be  accompanied  by  device  thermal  design,  i.e.,  the  choice 
of  dimensions,  materials,  and  processing  techniques  that  enhance 
heat  conduction  within  a  few  micrometers  of  the  device.  This 
work  provides  a  basis  for  the  thermal  design  of  SOI  FETs. 

This  work  shows  that  the  steady-state  channel-temperature  rise 
in  SOI  FETs  due  to  Joule  heating  is  significant.  Electrical-prop¬ 
erty  measurements  performed  on  devices  in  the  steady  state  are 
affected  by  this  temperature  rise  due  to  the  strong  dependence  of 
the  channel  mobility  on  temperature.  The  steady-state  data  may 
not  be  applicable  to  devices  in  an  integrated  circuit,  where  the 
channel-heating  is  lime  dependent.  More  work  is  needed  to  de¬ 
termine  time-dependent  temperature  fields  in  integrated  SOI  cir¬ 
cuits. 
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